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A CORRELATION ANALYSIS OF 
BOUNDED SOUND FIELDS 
by 

CHARLES JOSEPH GLASS 



Submitted to the Department of Electrical Engineering on 
August 23, 1954 in partial fulfillment of the requirements 
for the degree of Master of Science . 

ABSTRACT 

This thesis is an extension of the more comprehensive 
study made by Kenneth W. Goff (5) of the applicability of 
correlation techniques to the field of acoustical measure- 
ments . 



The ability of correlation techniques to separate an 
acoustical signal at a point into components according to 
(1) the source (2) the transit time from the source to the 
point in question and (3) the frequency; was used to 
measure the transient response of a reverberant room. The 
results of this correlation analysis were compared with 
pulse measurement data taken under similar experimental 
conditions. The results of the two methods were found to 
agree . 

A technique for measuring the directive properties of 
the wave fronts associated with sound fields was devised 
and tested. It wa 3 shown that the property measured cor- 
responded with' the definition of "diffusion.” Results of 
diffusion measurements taken in a reverberant room with 
various portions of the wall covered with a highly 
absorbent Fiberglas curtain were presented together with 
corresponding transient response data. 

The above mentioned experimental measurements were 
paralleled by an explanation of the theory and a sample 
of the mathematics involved in the correlation analysis 
of bounded sound fields. 
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I INTRODUCTION 



1.1 III TORY 

Tho word "correlation" inplioo tho act of cochinc 
rolatlonnhlpc botvocn two dlocroot quantitioc or functiono 
of tlno or noro proporly, two onsombloo of events. Tho 
concept ic olonontary . For that roan on, It la inpooclble 
to pin-point oono oxanplc and cito it aG boiriG tho flrot 
lnotanco of tho application of correlation. It io ouf- 
flolont to cay that tho procoos in ono which oiioo itn 
development to nan'o innato curloDity coupled with hio 
ability to roaoon. 

Historically, correlation 1c a vory important pro- 
cess. Mien nan boce.no diocatinfiod with tho noro 
rocordinc of ovente and first attempted to oxplain then 
in tome of a preceding carlo a of Incldcnta, he v;ac ro- 
cortinc to a fom of correlation. 

Correlation, liho nany other prococcoc involving 
basic concoptc, londa itcolf to concico and explicit 
oxprocclcn by mathematical oynbolc . The translating of 
tho procooc froa uorda to nathonatical cynbolo noro or 
loco parailolc tho ©volution of tlio procooc from an art 
into a ccionco. 

Tho natlionatical oxproooicaa of tlio correlation pro- 
cooc wac a natural outgrowth of tho dovolopoont of tho 
field of ctatlctical methods. Before 160C, no 



mathematical conceptions of probability voro rococnlred. 
Guilford, In hlo book on psychometric methods'*’ cays, 
"Gamblers had opoculatod nuch c once mine canes of chanco 
when It cano tlno to conoldor their Iosgoo and edns, 
particularly tliolr loonoc. They ovon attempted to intor- 
ost nathcmaticlcno In tholr problems, thouch with mall 
ouccoss. ... Tho covonth century saw tho boclnnlnc of 
oerlous lntoroct In the nathornticn of chanco. ... 

"Bernoulli (1634-1705) publichod tho flrct book to 
bo ontlroly devoted to tho subject. Dc!*lolvro (1667-1754) 
aay be crcditod with tho discovory of the normal diotrlbu 
tlon curvo at about 1733. Fran that tire on, lntoroct wa 
aroucod amonc aotronenoro as troll ao mathematicians . By 
1012 Laplace (1749-1027) had written that 1c cone Ido rod 
tho c^atect oinclo work on probability. In It ho qq.v o 
proof of tho method of leant squares. 

"It was Gauen (1777-1054) who denonctratod tire Groat 
practical value of tho normal curve, ohairine how it ap- 
plied to tho distribution of noacurcnontc and to errors 
node In scientific observations, it was ho who dovlsod 
tho fundamental noth ode of ccnputation of reans , probable 
orrorn, and tho ll!:o. ... 

"The application of tho normal curvo and elementary 
statistical methods to blolocical and social data must bo 

Roforenco nrabors refer to tlio Bibliocraphy. 
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attributed first to Ouctolct (1796-1079), royal astronomer 
to tlio bine of DolGiva. ... ir Francis Gal ton (1822- 
1911) in working on tho problems of hunan heredity found 
tint tho noiml curvo and its oinplor applications were 
inadequate. Ho lnvontod a nunbor of additional statistical 
tools, anonc thorn tho nothod of correlation." 

In 1CC8, olovon years after Galtcn’s firot publica- 
tion of tho concopt developed in torno of tlio ideas of 
rocro scion linos, tho tom "correlation" first appeared in 
print. Pearson, E&go worth, end Golden furt-ior dovolopod 
tho mathematics of correlation to a point where it was 
Gonorally adopted by tlio statisticians in the fields of 
economics, socioloGy, bloloGy, aotronooy, notoorolocy end 
a fow othor physical sciences. Hold, in Me booh on 

O 

statistics** reviews tho various applications of tho tlioory 
of correlation. An an addendum to his treatment of tho 
subject. Hold compiled a blbliocraphy which is excellent 
in its scope of tho subject of correlation. 

Perhaps tlio most important advance in tho fiold of 
statistics has boon tho mathematical exposition of the 
relationships botwcon tho fiold of complication ciiGlnoor- 
inc and statistical treatment of time series. Tlio rlcid 
mathematical proofs of Herbert Wionor^ not only put the 
procooc of correlation on a firm mathematical foundation, 
but also tended to ac colors. to further advancomonts in the 
fiold of communication onGinoerlnG . 
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1.2 TUI, 3KRIE5 AMD COXtfXATION 

loner do fined time series as "Ooqucncco, die ere to 
or continuous, of quantitative data assigned to specific 
moments In tlno and ctudiod with respect to tho statistics 
of thoir distribution In tlno. They nay bo simple, in 
uhleh caoo tiioy consist of a single nunorlcally Given 
observation at oacli nonont of tho diocreto or contlnuouo 
baco sequence ; or nultiplo, in which caoo thoy conoiot of 
a number of no para to quantitioo tabulated accord Inc to a 
tlmo common to all." 

Tho application of statistical nothodo to comnunica- 
4 

tion engineering roqulroo a model of an explicit function 
of tlno those proportloo con bo doocribod in tozmo of ito: 

1. mean value 

2. mo an- square valuo 

3. pernor spectrum 

4. correlation function 

Tho function of tlmo, be side a bo ins capable of description 
in tome of a cufficiontly complicated cot of probability 
donsitioo, must alco bo a stationary process ; that is, one 
in which none of tho probability dlctributiain vhlch 
doncrlbo tho pro coos change with tino. 

It is necessary at this point to define the statis- 
tical quantities associated with tho measurement of random 
functions of time such as tho thermal voltage of vacuum 
tubes. 



Tho nean value: If f* (V) dV Id the probability 

that V(t) Hoc bottToon V and V + dV thon tho 
ncan value of a nolo© lo 



whore the notatic n <( y noano the ctatlstlcal average, 
as conputcd fron the probability donedtioo 'Which dcncrlb© 
tho onnenble. 

Tho time -avc raced noon value of a noico voltaco 
V(t) io 



A random function nay bo ascxnod ■ to vanich outcldo 
of tho tloo interval -T/2 ^-t — T/2 In which caco 





under tho rootriction that 




ISio nean-oquaro value ic 







Tho noan-oauaro tine -a vo raced value io 





V(t) = 




o 
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where A(f) lo n cooplo:: function and oqualc tho 
voltage ooectrun of V(t) . ith tlio aid of 
Parceval * c tho ore n : 



| V s (t)dt 

-cO 



-Vz 

V°(t)dt * \lMt)l 

-Vz 



df 



wo nay obtain an oxproosion for tho total (finite) powor 
of a rand on function: 



lin _L 
tv*® T 



L 



Vz 



V a (t)dt = v*(tT = 

-T/z. 




df 



whoro r ./(f) , tho power spectrum, io defined as 



lin S /A(f )/ 8 

The auto- correlation function (p,X^) of a randan 
function lo 



= v(t) v(t-r) * <V(t) v (t -r)} 



similarly, the crooc -correlation function of two 
random functions io 

QtV = V x (t) V* (t- r ) * <(v x (t) Vg (t-r)^ 

The linlt botuoon tlio correlation function and tho 
other statistical quantities used to doocrlbo random func- 
tion has its basis in ./loner’s thooron: tho auto correla - 

tion function in tho cosine Fourier transform of tho 
intensity s-x>ctrun « 

Tho General relationships between the function of 
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tlno, tho spectral doncity, tho poiror cpectrun, end tho 
correlation function aro ah own in a schematic lo vised by 
Kbrmoth . Goff”*, Pig. 1.1. 

Tho results obtained by tho correlation analysis of 
acoustical problems aro similar in sono rospoctc to rocultc 
obtained by both otcady otato methods cvnd by pulso notholo. 
Fren a rudinontary imoulodgo of tho nature of spooch and 
nuaical coundc, ano goto tho intuitive fcollnc that, ainco 
tho sounds normally cncountorod in acoustical problona 
rocomblo pul nod \ravo trains, streosinG tho pulso-libo 
an jo etc of correlation no thole would yield slgnif leant 
rocultc, However, cuch a supposition io not cufficlont 
to justify tho nccloot of tho steady otato aspects of tho 
correlation process . 15 . io ability to apply a do ere o of 

frequency filtcrinG in tho cource of computing tlio cor- 
relation function has ccvezml advantages . Paramount among 
tlio co covoral ad van taco a ic tlio fact that much of tho 
present body of knowledge on acoustical problems io based 
an toots node with stoady ctato, frequency filt-orod 
noises. Tills fact end tho consequent advent aeon G^-i^od 
by tho ability to compare results with those obtained by 
noro classical methods make it desirable to conpronico 
tho time filtorinc properties of correlation methods. 
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Pig. 1.1 Interrelation between two time functions and their amplitude 
spectra, cross power spectrum and crosscorrelation function. 
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i,3 ron niCTsan vs nsnracx fiIttriico 

I !ony of tho &':nlicatlcaia of correlation tochniqucs 
in the fiold of acoustics do pond upcm tho ability of th© 
nothod to separate signals into ccnpomonts according to 
thoir transit tlno from a source. 

For ins ton co, suppose two micro phono a xioro placed in 
tho sound fiold of a point source radiating randen noiao 
into froo space. If tho nicrophono outputs wore nlrrod, 
and this conbinod voltage wro crooc-corrolatod with tiio 
loud spooler input volt ago » tho follotfing of foots would b© 
noticed: At a do lay tine corresponding to tho transit 

tlno frea tho loudcposicor to the microphcaia closest to 
tiio sourco, tho correlation function would exhibit a pools 
value. At a do lay tlno corresponding to tho transit tino 
fren tho loud spooler to tho cocond nicropiicno, another 
poah would appear in tiio plot of correlation function vs 
tlno delay. 

If tho nicrophono s were placed too close together, 
tho rosultont plot of correlation function vs tino delay 
would Ido distorted. Tho smallest difference in transit 
tim which producoG an undistortod nain peak in tiio cor- 
relation fine ti cn corresponds to tho ro solving poorer of 
the system . Tills nininun transit tino my be of tiio ordor 
of ono period of tho geometric no an frequency of tho band- 
width of the random nolso used; and the corresponding 
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rooolvlng po’.ier of tho system lo equal to tho product of 
this minimum transit timo and tho cpood of sound In tho 
air. 

Since tho auto-corre lotion function and tho po^’or 
apoctrun of a stationary randan tino function arc Fourior 
transforn pairs, thoy exhibit certain characteristic in- 
vorco spreading proportios . For example , a froquoncy 
spoctrum, F(o) , flat fron 0 to - rt/r lias as its 
transform, a function of tino , f(t) , having tho general 
shape of 1 sin t which has one predominant peak and 

V ZL t 

r 

axis crossings at t T , t 21, - 4 T, ... , Fig. 1.2 

(a). As tho bandwidth of F(o) Increases, tho axis 
crossings of f (t) occur at xjroGrocslvoly smaller values 
of time. Tho limiting ease is an F(oj) equal to 1. Tho 
Fourior transform of such a spectrum is the so called 
"impulse" or "de/fa function" which is a function of time 
having an infinitely Gmail duration, an infinite ampli- 
tude and an area of 1 under tho curve ; Fig. 1.2(b). Con- 
versely, as tho bandwidth decreases, tho peaks in tho 

function sin x tend to equalize. In the limiting 
x 

case, F(a) is composed of one froquoncy (- o©) , 

and tho inverse Fourior transform f(t) is 

cosino (o©) t ; Fig, 1.2(c). 
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Figure 1.2(a) Fourier transform pair for finite band- 
width noise. 



f(t) 



F(oi) 



1 

1 

L 

o 



o 




Figure 1.2(b) Fourier transform pair for infinitely A 

wide bandwidth noise. 




Figure 1.2(c) Fourier transform pair for a pure tone. 
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From tho above discus clan and llluctratlon It is 
obvious tliat tho proroquiciteo for high resolution In 
tine aro diametrically opposed to tho requirements for 
high frequency rooolution. It Id equally cloar that noot 
acoustical noacurononto undo by corrolatian no thod a ro- 
quiro a comproal do botuoen tho frequoncy and tine 
filtering proportios of tho process. 

For instance, in measuring transmission loss of 
panels , one generally do sire a to laiov; ho u this property 
varies i/lih frequency. Efrus tho bandwidth of tho tost 
signal must bo cut to a dogreo commensurate with the fre- 
quency resolution desired. Eh is lncreaso in frequency 
roooluticn must not bo so groat an to croato confusion 
botv/oen tho ouccosolvo peaks caused by tho signal trans- 
mitted directly through tlio panel and those peaks caused 
by signals arriving at tho tost point by dlfforont paths. 
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II INSTRUMENTAL I ON 



2.1 IXTROWCTIGK 

Iho eubjoct of instrumentation must includo con- 
sideration of tho methodo of (a) g& tho ring data, (b) 
reducing data, (c) presenting and recording results. 

Equipment for Catherine data nay lncludo ouch itoms 
as oloctranic noise voltage Generators, loudspeakers, 
microphones, and in a caoc ohero the reducing equipment 
io not portable, recording instruments . A reducing system 
io gonorally conpocod of a number of separate oloctronic 
unite oach of which perforins a opocific process. Division 
of the reducing procoeo into a sot of sub-prococooo Ionia 
verr, utility to the eye ton and also facilitateo testing of 
the components. The end product of tho reduction process 
nay bo presented in various ways. If the reduction cysten 
io oloctronic, tlio rooulta may be presented by means of 
motor readings or oscilloscope traces. 

If a reducing system is already In existence, then the 
equipment for gathering data and presenting results must 
bo soloctod in such a way that the best uoo will bo made 
of the features built into the reducing system. 

2.2 THE ANALOG CORRELATOR 

As implied by the oouati on defining tlio croso- 
corrolatlon function 



21 ( r) s lim 1 

T-»* T 
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a device capablo of commuting this function nust porforn 
tho operations of (a) delay, (b) multiplication, and (c) 
integration. If tho functions of tlno vero oloctrical 
voltages ouch as thooo produced by a microphone In a 
sound field, the application of those voltages dlroctly 
to the input terminals of an oloctrical analog computer 
dosigned to perform tho required oporatlon would yiold 
tho correlation function dlroctly as an output voltage 
which could bo measured end recorded by properly nolo c tod 
instruments . 

Tho block diagram of ouch an analog ccmputor designed 
and developed by Kenneth V/. Goff ^ is shown in Fig. 2.1 . 

A x^iotograph of tho correlator io shown in Fig. 2.2. For 
tho purposes of analysis, tho computer Id composed of 
throo ooparato units 

1. tine delay system 

2. multiplier 

3. integrator 

Each of the so units nay bo oporatod independently of tho 
remaining components. 

2.2 (a) THE TIIIE DELAY 

The time delay mo eh onion employed in tho computer 
being canalderod consists of an olc-ctro-mechanical dovico, 
Fig. 2.3. A* shown, in tho block diagram, t!io two signal 
volt agon being usod are fod into separate channels each 
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Fig. 2.1 Block diagram of the c orrelator. 
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'Fig. 2.1 Block diagram of the correlator. 
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Fig, 2.2 



Ifcatograph of analog oorrolatim compute* 
cyctoa. 
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having its oim varlablo Gain ompllfior, magnetic rocord- 
lnG hood, magnetic recording track, playback ho ad, ora go 
hoed, and playback amplifier. Tho record head of cimnnol 
one in fixed, whoroao tho record hood of channol two is 
movable. If tho transit timeo from tho rocord to tho 
playback ho ado for a point on tho ourfaco of tho roc online 
drum are equal for both channels, then no relative tlmo 
delay ’id 11 bo introduced between tho two input voltages. 

If tho transit timo between rocord and playback hoado lo 
groator on channel two than on channol one, then a rela- 
tive timo delay equal to the difference in transit timoo 
is introduced. For identical input voltages to oach 
channel, tho output of channel one would bo f (t) and tho 
output of channel two would bo f (t-T) where T io the 
difference in transit time. 

This relative time delay betwoon chamois ono and 
two can bo varied continuously from -15 milliseconds to 
190 milliseconds by moans of gear train and chain drive 
system powered by a low speed synchronous motor. 

Tho frequency response of tho system io flat to 
within - 2 db ovor a frequency range from 100 epa to 
10 keps. 

2.2 (b) THE MULTIPLIER 

The multiplier componont of tho correlator is a 
quarto r-difforonco-cquarlns do vice which reduces tho 



procooc of multiplication to that represented "by the 
relation 
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fi fa = | (fx + fa)° - (fx - fa) a 

The process of producing a voltage proportional to 
the square of (fx 4- fa) and (fx - fa) ie per- 
formed by a squaring circuit consisting of tv/o 6B8 
pentodoo operated with grids in push-pull. For opera- 
tion in the region whore tho trano conductance and grid 
voltage ax’© lineally related, the plate current for ono 
tube can bo approximated by 

i, = a + ho^ + c© ® + de * 

d 6 6 6 

For the two tubes operating in parallel with a common 
plato load as shown in Fig. 2.4, the current through tho 
common load is tho sum of tho plate currents 

*l = V + ip" = aa + a™! 0 

The output voltage 

6 o= E bb- 1 L E I,= at V 

If two squaring circuits are operated in push-pull 
and e^ is made equal to (f x + fo) in. ono squaring 

and equal to (fx - fa) in tho other squaring circuit. 




£ d ^iiripiifieci Sv hema*i( of the q lann^ ircuit. 
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then the output voltage will bo proportional to the 
product fi f® 

Variability in tho performance of the correlator 
generally results from unbalanced conditions in tho 
multiplier ciroulta. When tho multiplier in opera tins 
correctly, tho output voltage should bo aero when there 
ia no input and uhen there la only ono input. Tho so 
balanced condltiono correspond to tho requirements that 
OxO ss 0 and tel = 0 = 1x0. 

Tho balance for 0x0 = 0 is made by adjusting tho dc 
bias of the four 6138 tubes for equal quiescent plato cur- 
rents whereas tho balance for tel = 1x0 = 0 ia made by 
adjusting the ac levels applied to each 6B8 tubo. 

Much attention was focused upon tho design, require - 
mento no co scary to reduco dc drift to a degree that would 
permit continuous operation for periods up to two to 
throe hours without adjustment* . 

2.2 (c) THE INTEGRATOR 

Integration in tho computer is accomplished by moans 
of either a low pass filtering system or a stopping 
integrator. For acoustical measurements, a simple RC 
low pass filtering circuit is satisfactory. Hie corre- 
lator has several such circuits with RC timos of 0.5, 

1.0, 2.0, 4, 8, and 16 seconds. Selection of tho proper 
RC time is another important factor governing optimum use 
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of tbo equipment . 

The dc ov.trrat of tho RC circuit is amplified pnd 
then placod across tho input of a dc-ac chopr>er \/hlch 
converts tho dc voltage into a proportional 4 00 cor 
signal . Tills 4 00 cps signal is tlion amplified again and 
placod across the output terminals of the integrator. 

Tho 400 cps signal may be measured directly with a 
vacuum tube voltmeter or an oscilloscope, or it nay bo 
rectified and recorded by moans of a dc graphic ammeter. 

A linear plot such as is obtained on a graphic ammeter in 
extremely usoful since it represents tho shape of the 
computed correlation function. Information regarding tho 
relative amplitude in decibels of successive peaks of tho 
correlation curve i3 desired. Such information nay bo 
obtained by moans of a logarithmic graphic lovol recorder 
which ro cords the correlation curves directly in db. 

Yfoon usod together, the linear and logarithmic plots give 
a continuous, detailed, visual record of tho shape and 
magnitude of the correlation curves as a function of time 
delay. 

All of the information presented in the foregoing 
discussion of tho analog correlation computer has boon 

g 

taken from the doctoral thesis of Kenneth V. Goff . 

Only those points which are considered necessary for a 
basic understanding of the computing procoos are pre- 
sented here. For a more detailed treatment of the 



docign roqulromcnt and tho performance characteristics 
of tho computer, tho reader is roforrod to either rofor- 
cnco 5 or to tho recently couple tod, unpublished namnl 
"Instruction Ilanual for Analog Correlator" by Kenneth W. 
Goff, ..ugust 1954, M • I . T . Acoustics Laboratory. 

2.3 DATA GATHERING ACCS33CRIS8 

Tho accessories selected for gathering data aro as 
follows : 

1. The noise cource consisting of an electronic 
noiso voltage generator, a power amplifier, and 
an Altec 728B loudspeaker in a closed box baffle . 

2. Tho microphones used were Altec Lancing typo 
21 BR-150 wi th 157 -A bases and P518A Altec 
Lansing power cuppllen. Tho response curves for 
these microphone 3 are flat to within - 2 db for 
a frequency range of from 100 cps to 10 keps. 

The voltage amplifiers uced with these micro- 
phonos wore wide band amplifiers. 

3. Since tho reducing equipment was not portable, 
provision for recording data v/a* required. An 
Ampox model 350-2 twin track magnetic tapo 
recorder (Fig. 2.5) ’was uced for this purpose . 

The particular apparatus used had a frequency 
response flat to within - 2 db ovor a range of 
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FlS* 2*5 ^totograjfa of mpox nagnetic tone rocorGor 
vdth a tap® loop. 
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150 cps to 15 kcpo. Ella oqulpaent v:ac 
portablo and yot capable of reproducing a 
relatlvoly dlctortlonloaa roplica of tho In- 
put voltage. 

2.4 ACCESSORIES FOR PRESENTING RESULT^ 

The consolation function as computed by tho analog 
correlator Is pier-ented in tho fora of a 400 epe voltage 
and a roctlflod dc voltago. It is therefore possible to 
measure and record tho correlation function by neons of a 
variety of electrical measuring instruments. 

Tho accessories ucod for no a curing voltages arc: 

1. Ilodel 300 Ballantlno vacuum tube voltmeter. 

2, Typo 208 Dumont Cathode Ray Oscilloscope. 

Remanent records of the output voltage x.’ore made 
by noans of: 

1. Model A'./ linear scale Esterlino Anguc Graphic 
Lilllamnotcr. 

2. Bruol and Kjaor log scale Graphic Level Recorder. 
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III. TRANSIENT RES PON Cl OF ROCIIS 
3.1 INTRODUCTION 

Ilony studios ha vo boon nado on tho subject of sound 
fields in rooms. Generally those studios havo boon based 
on ono of two distinct approaches to tho problem of room 
acoustics : 

a) a study of tho effect of rocans 1 boundaries on tho 
sound field associated with steady ntato pure tone 
signals. 

b) a study of tho effects of rooms' boundaries on tho 
sound fiold associated with transient signals. 

Thoroforo, in attempting to obtain a complete quanti- 
tative assessment of the offocts of boundaries on a Bound 
field, ono mi$it: 

a) excite the room with a non-directi anal s toady state 
sinusoidal signal and map tho amplitude and phase 
relations of tho sound proosuro for 

1. all frequencies 

2. all possible observer positions in tho room 

3. all possible source positions in tho room 

b) oxcito tho room with a pulso of oxtremoly short 
duration and measure tho relative amplitude and 
diroetian of prorogation of tho associated wave 
fronts for 
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1. all tirao 

2. all obsorvor positions 

3. all source positions. 

If ono has all of tho data associated with oithor 
ono of tho abovo mcntlonod oxporlnontal approaches, ho 
can mathematically deduce all of tho Information for tho 
other method. For example, ono can synthesise tho pulse 
used In tormo of a Fourlor serios of oinuooidals. Then 
for oach frequency component ono can pick tho rolativo 
amplitude and phase angle of the sound preosuro at any 
point in tho room. By weighting those amplitude and phase 
measurements in accordance vrith the Fourlor analysis of 
tho pulse and then adding the effects for each froa.uoncy, 
tho complete response of tho sound field to a transient 
pul so can be deduced. 

On tho other hand ono can visualize a model of a 
room excited by a pulso. In ouch a model, tho walls of 
tho room would bo eliminated and thoir of feet on tho sound 
field would bo simulated by an infinite multidimensional 
array of imago sources, oach of which is located in a cell 
of tho samo si so and shapo as tho actual room. Each imago 
source is positioned so that tho transit tine of tho sig- 
nal from tho imago source to tho actual observer’s 
position equals tho transit timo for tho signal arriving 
at the observer’s position by a particular sequence of 
reflections from tho room's boundaries. 
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With thlo model, ono co.n iw/ specify that oach non- 
dlroctional source is In phase and is radiating a froc 
travel Ins sinusoidal signal . By adding up tho contrlbu- 
tian from each imago oourco, ono can nap tho amplltudo 
and phaoo relations of a room excited by a o toady state 
source. 

The above discussion assumes that the boundaries of 
the roca are inflnitoly hard go that no phaso ohlft or 
do ere a co in amplltudo occurs when a signal lo rofloctod 
from the walla, Thlo assumption slmplifioo tho problem, 
but tho ro a ul ting calculations do not conform with 
reality. It is nocoasary thoroforo to consider those 
phase shift n and attonuaticna occurring at tho boundaries 
of tho cound field. In the transient treatment of aound 
fields ono may account for the boundary effects by as- 
signing to each bounding surface an impulse rosponse 
H(o) , In ctudylng tlio steady state behavior of bounded 
cound field a an absorption coefficient and a phase shift 
are assigned to oach wall. 

Any ouch complete analysis of a bounded aound field 
ifould be extremely tedious. A deeiro to bo practical 
would incline ono to restrict tho analysis to a degree 
commonsurato with the amount of detail no cocoary for tho 
proper utilisation of tho room. In moot rooms we are in- 
terested in producing an acoustical environment conducive 
to good audl o communications. Therefore wo may generally 



29 



ro strict our consideration to tho audible frequency 
ranG© . Audio c orrauni c at i ann will generally fall into 
one of threo catog orien. 

1. two way vocal c ocrauni cations botrecn individuals 

2. vocal communi cat ions between one speator and an 
audience 

3. instrumental cocmmi cations between musicians and 
an audience 

Traditionally, each type of communication mentioned 
above has associated with it a certain physical arrangement 
of source and observer. If we accept those traditional 
arrangononts as bo ins characteristic of the type of com- 
munication to be used in tho room then wo may further 
limit our analysis of a bounded sound fiold to those 
source and observer positions of practical into rest. 

For practical purposes we may further restrict our 
analysis by considering only those types of sounds char- 
acteristic of musical or vocal sources. Such simplifica- 
tions further restrict both tho frequency composition of 
the ton tins signal s and pulse duration and pulse inter- 
vals for transient signal to sting. 

Up to this point we have been discussing tho two 
general approaches to tho study of room acoustics end tho 
various simplifications that may be employed in this 
study. Those two approaches have been quite thoroughly 
explored and tho resulting body of data coupled with 
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subjective evaluations of tlio acoustical environment In 
rooms has lod to tho conclusion that variations In both 
tho steady state characteristics and tho tranoiont char- 
acteristics of sound flolds account for difforonco3 in 
tho acoustical quality of a room, 

'/.C. Cabin©'’ recognized tho import anc© of the 
transient behavior of sound fiolds and attempted to re- 
late it to a subjective evaluation of tho quality of 
roomo. As a noasuro of transient response, Sabino 
defined a * ro vo rbo rati an tine” the length of time for 
tho no an square ores sure of a suitably chosen distribu - 
tion of sound v;avos to diminish to ono millionth of its 
original intonsltv . Rove rbo rati on time as dofinod by 
Oablno Is a me a euro of tho rato of docay of the sound 
pressures caused by a steady state sinusoidal signal 
which is suddenly discontinued. Various criteria based 
on rovorbo raticn tlmo vs tho si so of a room which is to 
bo used for a specific application, and rovorboration 
tine vs frequency have boon devised 1 . At present, 
design curves of this type offer tho most satisfactory 
approach to tho construction of rooms having good 
acoustics and tho treatment of rooms having poor acous- 
tical characteristics. However, it is recognised that 
acoustical design crltoria other than rovorboration tine 
must bo considered in tho correct design of rooms. For 
instance, ono research group^ suggested that tho first 
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20 db of docay of a cotmd in a rood io of primary im- 

portanco in differentia tine bcf.ro on two rooms which have 

approximately oqual over-all rovorboration tlnos. 

o 

C.A.. Hasan and J. Uolv * , in the course of their 
studios of tho acoustics of theaters by means of pulses 
or '’tone bursts , n Intimated that tho tine and anplitudo 
distribution of reflected tone bursts might servo ac an 
aid in evaluating tho acoustical quality of rooms. 

In tho introduction to "Pulse Statistics \nalysio of 
Roon Acoustics" tho authors state, ”... Thoso facts indi- 
cate that . . . tho first few reflections aro primarily 
responsible for certain important features of the acous- 
tical character of rooms ... It should theroforo bo 
worthwhile to study this ‘short torn* transient ro opens© 
by a method of images in which oil of tlio wave properties 
of tho image sources can bo considered (i.o. whore tho 
assumption of an Incoherent source is not made). Further, 
if on imago array satisfying the boundary conditions can 
bo found, one should be able to treat this array statis- 
tically and thus obtain tho long torn average transient 
response as well ...” 

As indicated by the above quotation, sene method of 
measuring tho transient response of rooms io required for 
on evaluation of the effects of boundaries on a sound fiold. 
However, ccno estimate of tho steady state behavior of 
rooms is equally desirable. One experimental nothod of 
measuring the transient response of roams lias been 
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GUGCOotod .iiich oorao to combine some of tho pro; rtios 
of both steady t to raothoae and pule© methods’ . Thlo 
method employs correlation techniques. 



3.2 7 in COP .EL ‘TI N YSALYOiq 07 TTAK'ITNT ^SIGH'S: 



How that we havo suggested cooo potentially useful 
mothodo for making acoustical measurements, it would be 
wise to explore briefly the mathematics of thio mothod. 

It can be oho'i/n directly that tho correct applica- 
tion of correlation measuring tochniquos to bounded sound 
fields will yield the impulse response of that bounded 
sound field . 

Hof erring to Fie. 3.1, w© as sumo that the noic© 
source is do live ring a signal having the complex spectrum 
F s (to) to the two channels of the correlator. If tiro net- 
works havo transfer functions Hi (2) and Ha (3) respectively, 
tlren tho output of tire networks aro : 

F a (o) =s H x (Ja>) F Q (w) (3.2-1) 

F b (cd) ss H«(Jc d) F 0 (o) (3.2-2) 



Tho cross power spectrum for Fx (to) and F 0 (o) is 



$ 01 (cs) = lin rr 
T 




F s (o) Ha(-Jai) F a (d)] 



= $ GG (a) Hi(Jo) Ha (-jo) (3.2-3) 



CORRELATOR 
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■where 



= 



lln tt 



£ 



(o) ) 



(3.2-4) 



la tho power density spectrin of tho nois© source. 

Now tho correlation lo tho Fourior coa5.no transform 
of tho power density spectrum 



<P 
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(M = 




(©) Hi (Jo) H«(-jo) O' 1 do 

(3.2-5) 



Frea Fig. 3.1 t/o can noe that frequency variation in 
the source can bo considered as common to ran to both chan- 
nels one and two. This factor then can bo considered as 
a constant, i.o. 

<P 21 ( T ) • $ oa f h a^ q ) Ha(-Jca) e ^ a> dm 

( 3 . 2 - 6 ) 



If we now consider the opoclal case ’.die re H s = 1 
and wo wish to auto-corrolate tho output and Input of tho 
system having a transfer function of Hi (jo) , tho croas- 
corrclaticn function lo 



21 



( r) = § 



L 



Hi (jo) 



jar 



do S5 



art $ aa hi (t) 



(3.2-7) 



Thuo, the cross-correlation function lo exactly 
proportional to tho Impulse response of tho system having 



35 



tho transfer function Hi (Jo). If tho channol having tho 
transfer function Hi (Jo) la a roca, then the correlation 
function 5.0 proportional to the impulse roaponco of the 
roan. 

'o havo talked a little about the nature of tho cor- 
relation function; lot us now concider lta shape for a 
correlation analysis of tho sound field of a room. For 
this purpose, wo shall first define a fow terns . 

Using Laplace Transformation wo mdo tho definition 

Pi (o) = 1 f fi (t) e ‘■ 1cat dt (3.2-0) 



iSionco 

fa. (t) ts I F x (o) O Jcot do 

The ITultipllcotion Theorem states: 

2tt j^Fx(o) Fe(o) G cr do = 

f fi(t) f# (t-M dt (3.2-9) 

Now dividing each cido by 2T where T a tlmo, 
and then taking the limit as T 



lim 2tt 

T->*° 



j Fi(o) FoTa) e * ot do = 



lin 



.1 

2T 




fi(t) fo(t-T-) dt ^ <^ 21 (D 

( 3 . 2 - 10 ) 
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In tho application of correlation techniques to tho 
trc.no lent ronpon.ee of rooms fi(t) in tho voltage produced 
at tho loudcpoabor terminals by a vrido load noise-voltage 
generator and f Q (t) io tho voltage output of a microphone 
placod at sono observation plaoo In tho sound field. 

Tho volt ac© produced by the microphone io fe(t). It 
c anoints of tho cun of fi(t) nodiflod by tho transfer char- 
actoriotic II of tho loudspeaker and microphone, plus a 
numaticn of reflections ©ach of which oquals f x (t) 
dolayod by a tino equal to tho path length of transmission, 
1, divided by tho opoed of sound Cl; and modified by H and 
a transfer charactoriotic B aocociatod v,lth its transmis- 
sion path. 

The frequency domain representation of such a voltage 

is 

je> iJo io lx/c 

F« JB Fx II B 0 o * + F x II Bx O + 

jo 1*3/0 

FxIIBa e + (3.2-11) 

The B’o oro ccmplon coofficionto which account for tho 
phase shifts at tho reflecting surfaces, tho attenuation in 
amplitude associated with Invoroo square dizsunlticn of 
radiated sound Intensity, and tho attenuation in amplitude 
occurring at tho ro fleeting surface . 

B * IbI e Jo>P 
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thorcforo 

Fa =: 



nO 




n i 



-Jo "K/o 

F D H o ^ ^ 



(7.2-12) 



Ci J (^ - ol^c) 

= V ^ l°nl ° 

>1'- i 



If e now nobe the definition: 



c rs poorer density spectrin = lln 



then It follows that 



F(o) 



= u a TEi 

w** / rr 



e 






rr 1 f (o)/ * 

T 

(3*2-13) 



(3.2-14) 



where 



F(to) = 



/f(o)/ 






low substituting in the nultlpllcrvtion theorem: 

<P v2 ( T) s lia 



“3 (^1 “^9 ♦ ot) 



I n 75 j 

If / rr 



^ LSi 

rr 



do 



a jVsiSj e * (ot -<* x ) (3.2-15) 

'..lion, we correlate F x and ?c using equation 3.2-15 
and 3.2-12 



*12 < T > = 



if ^ e 1 3 n 



i <P 



n 



- oyell 



do 
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<P 



12 



(M * 




0 J <p " + B<r ' 1 n / «^ - 
(3.3-10) 



Equation 3.2-16 then Givoo tho crooo -correlation 
cp 7" ) between the Input voltage to tho loudnpoa 1 :er 
end tho output voltage of tho nicro phono in to ran of tho 
power density spectrum (3) of a noise voltage Generator, 
tho combined trannfor function (H) of tho loudspeaker 
and nlcrophcno, and the complex trannfor function of tho 
transmission path. Chic lattor transfer function 
/ V 0 }P ' n dor ' crtlx,a em ^ x modification in tho 

sound pressure amplitude c poet ran caused by its traveling 
tho distance involving any number of reflect! ana. 

For tho idealised caso where: (a) nourco radiates 
noico of a flat cpoctrun from a x to oo and soro outaid© 
of this ape c trim; (b) tho transfer responses o of tho loud- 
speaker, microphones, and walls aro 1 LQ” } it liao been 
TO 

shown that tho crono- correlation curve for tho loud- 
opealzor and microphone voltages will have equal pc alia 
for cadi point where tho tin© dolay corresponds to the 
transit tine fro n tho source to tho microphone and that 
tho cross-correlation curve v&ll damp out to a email 
value be tiro on tho so points. Tho rapidity with which It 
damps out will do pend upon tho bandwidth (©a - oiJ and 
tho torn (o a + oi). 
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with the foregoing mathoaatica in mind we may no\/ 
visualise a room model that will bo useful for correla- 
tion analysis. Thic bo&©1 will be exactly the oamo ao 
tho nodel mentioned previously In section 3*1 in conmc- 
tlcn with puls© analysis of rooms. In this nodel the walla 
are replaced by an array of imago aourcoo each of ;:hich la 
located at a distance equal to the length of tho trancnlc- 
aian pa tli of tho reflected sound. Tho output of oach 
oourco nay bo nodi fled in accordance with the nunbor of 
roflccticna that occur in tho path of transmission from 
tho ocurco to tho obsorvor’s position. 

In their discussion of puloo statistics theory, 

Bolt, Beak, and vrootorvolt^ dovolopod tho mathaaatics of 
ouch a nodel. In their model oach imago oourco occupioo 
a go para to imago replica of tho phyoi cal room. Those 
imago roo os or colic aro designated by throe numbers 
(1> n, n) each of which can talas all integral values from 
minus infinity to plua infinity. Tho notation deolgnating 
tho imago source a is also an indicator of tho nunbor of 
rofloctlonn that tho signal undorgooo in traveling from 
tho source to tho observer, 

Vo shall now mate a few simple adaptations of thio 
spatial model to aid us In visualising tho correlation 
analysis of tho problem. All of the oourcoo shall bo 



Grouped in ac coni an co with tho number of reflections 
that occur in transit from tho courco to tho obcorvor. 
Thus the courcco uhooo sound io refloctod tvdco in trans- 
it aro docicnated ao oocond ordor inaQo courcoo. 

Tho nunbor of inaco courcoo in oach ordor dopondc 
upon tho room cooootry. For oxanplo, tho nodcl for a 
root angular roan can bo broken down into it a inneo coa- 
pononto follows t 



SOPRCSS HUnBIiP. 

Actual 1 

lot Ordor Inac© 6 

2nd Ordor Imago 18 

3rd Ordor Inaco 38 

etc. 



An outline of tho value c oncloeod by tho actual pluc 
tho oix flrct order inaco sources for a roctancular room 
io shown in FIs. 3*2. Tho nunbor of inaco courcoo of each 
ordor can bo computed from tho formula 

/V-/ 

nunbor s 2 + + 8 ^ k (3.3-1) 

where H lo 'tho ordor of tho reflection and k = (N-l) . 

For example , for tho fourth ordor array, there will bo 
2 + 16 + 8 (3*2+1) = 66 image sources. 

Other factors of importance in correlation cnalyoic 
aro tho distanco between inaco sources, and tho distance 
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Fig* 3 »2 Volume oncloned by actual and flrct order 
ini!i£j>0 o outgo □ for & rectangular roocio 
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from the imago nourcoo to tho observer's position. 

Thoco factors nay bo computed from fornulao 9 and 11 
of roforenco 8. 

The advantages gainod by reference to this spatial 
model of a room aro threefold? 

a) tho system gives a visual prooentatien of the 
problem o encauntorod in opaco or time filtering. 

b) the re cults of tho correlation analysis nay be 
compared with thOGo of other researchers who used 
pulse method 3 of analysis. 

c) the mathematics for pulco analysis dovolopod by 
tho use of this model may bo appliod directly to 
correlation analysis. 

3.3 EXPERIMENTAL MEASUREMENTS 

The experimental measurements of tho transient re- 
sponse of rooms wore taken by cross correlating between 
tho input voltage to a loudspeaker and the output of a 
microphono located in tho roam. A hard piaster room. 

Room 10-390-A, which had boon previously analyzed by 

O 

puloo methods was selected for tho preliminary tests. 
3oth tho oourco and microphone, as well as tho tost 
positions in tho room, woro selected to give results 
comparable to those obtained by Bolt, Doak, and 
Westcrrolt 8 . 
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Tho results for ono particular tout run aro shown 
in Fig. 3.3. Shown also aro tho comparable ro cults of 
a pul co analysis. In this m, tho source, a Model 175 
Jane 3 Lansing horn driving unit with a Model 1217-1290 
Jonas Lansing throat and lone assembly, was placed in 
the comor of tho room facing tho wall. Tho microphono, 
on Altoc 21-BR-150A vrac located 13 foot from tho acurco 
along tho short wall of tho room. 

Tho portinont dimensions aro: 

Roas dimensions : length Lx = 23 foot 

width Ly = 13.4- foot 

height Lz = 8 .4 feet 

x = 1.0 inches 
y = 1.5 inches 
z = 0.5 inches 
u = 1.0 inches 
v = 13 foot 
v = 1.0 inches 



Gourco Position: 



Mike Position: 



Tho tost signal was random noico gone rated by an 
electronic noise voltage generator and amplified by a 
Model 20-W2 McIntosh amplifier. A 1/3 octavo band 
analysis of the loudspeaker voltage and the nicrophono 
voltage is shown in Fig . 3.4. The loudspeaker voltage 
and the nicrophono voltage wore recorded simultaneously 
on separate tracks of a twin track tape recorder. In 
reducing tho data, a section of tho recording tape was 




Figure 3.3 



M1CM S^ED LIVtL RECOWOE" 
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cut and node into a tap© loop as shown In Pig. 2.5. 

Tho microphone voltage was filtered bo that its spectrum 
was that passed by a 1/3 octavo band filter with a Geo- 
metric neon froouoncy of 4000 cos. Tho filtered micro- 
phone voltage vrac fed into channel one of tho correlator 
tine dolay and tho tmf titered loudspeaker volt ago \;an 
placed across ch«annol two of tho correlator. An HC in- 
tegrator time of 1 second and a tine delay scanning rate 
of 1 nillioocond per minute were U9©d. Tho resulting 
linear plot of tho cross correlation function togothor 
with tho plot for the comparable pul no toot is chaxm in 
Fig, 3*3. In this figure, each group of peaks in the 
correlation function represents tho arrival of tho sound 
from on Image source in much tho sane way as the peaks 
in tho accompanying pulse analysis. An accompanying 
log plot of tho correlation was made. This plot is not 
presented because of the extreme physical length of the 
plotting tape. Ilowovor, Fig. 3*5 contains a plot of the 
peak values of the correlation curve as a function of C 
for Fig. 3.2. Plotted on tho same coordinates is a graph 
of tho peak pulse pressure level as a function of distance 
traveled. Tho points for this second curve arc talson from 
Fig. 4 (reference 8) . 



RELATIVE LEVELS IN DB 




FIG. 3.5 Plot of the transient roopcnco of 

Room 10-390A for pttlco end correlation 
analysis » 



3.4 DI5CUSSI0R 0? EXPERIMENTAL RESULT ° 

Tho data an pulso statistics analysis® shows the 
c pacing of pulses arriving at the microphone pooition 
when be til source and receiver aro In tho corner of tho 
room. Tho corrolatlon analysis of tho tranclont rospons© 
of tho room should Ghow comparable pooka In tho corrolatlon 
curve for tlno delays corresponding to tho no pulno spaclngo 
Tho correlation curvo, Fig. 3.3, shows discroto peeks at 
delay times of 12.7, 13.6, 21.2, 34.4, 37.2, end 47.5 allll 
seconds. Roforonco to Fi6. 14 of roforonce 8 choirs that 
those delay tines corro spend roughly with tho transit 
tines for tho various reflected pulses. 

Tho relative amplitudes of tho peaks of tho correla- 
tion give a quantitative noasuro of tho combined offocts 
of: 

1. Attenuation of oignal caused by invorso square 
dlninuation of radiated sound intensity. 

2. Attenuation of signal caused and decrease in 
coherence of radiation from tho various imago 
sources caused by 

a) Tho transfer function of one reflecting 
surface for the case of non-degonsrato 
discroto first order Imago sources. 

b) The combined transfer function of all the 
reflecting surfaces for higher order nen- 
dc genera to imago sources. 
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c) Decrease in coherence of tho radiation from 
groups of degenerate image sourcoo. 

Thus, tho plot of tho correlation function vo tin© dolay 
is a noasuro of tho short tora transient response of tho 
room for a particular sourco and receiver position. 

Using correlation analyoio, it is also poociblo to 
derive a value for rovorboration tiiao based on tho transient 
response of tho room. Tho rovorboration tine simply cor- 
ro spends to tho differonco in dolay timo botwen tho peak 
in tho curve corresponding to tho arrival of sound from 
source 0, 0, 0 and a peak in tho correlation curvo vrhich 
is ano -millionth of, or 60 db dovm from tho peak of tho 
direct sound. This tine nay bo found by oxtrapolating tho 
rate of decay indicated by tho short tern transient response. 
For tho example discussed in tho previous soction vhoro 
both sourco and rocoivor wore in tho comor, the correla- 
tion curve decayed 20 db in tho interval T = 10 nl ni- 
ce c and c to^= 67 nillioeccnds. The rovorboration tin© 
based on this rato of decay is 3.14 seconds. This value, 
then, is a noasuro of tho rovorboration tine for a 1/3 
octavo bond signal having a geometric moon frequency 
of 4cco cpo. 

In lino with tho experimental pulse analysis of 
reference 8, additional transient response data corres- 
ponding to tho data of Figs. 7 through 16 has boon col- 
lected and reduced by means of tho analog correlator. 
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Thin additional data uith explanatory no ton io preoontod 
in Figs. 3.6, 3.7» 3.8, and 5.9. For tho convenience of 
tho reader, a copy of reference G is included in tho 
Appendix. 




Figure 3.6 Top 
Figure 3.7 Bottom 
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Figure 3.8 Top 
Figure 3»9 Bottom 
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IV. MEASUREMENTS 0? DIFFUSION 
4.1 INTRODUCTION 

In the literature on room acoustics, much att rtlcti 
hao been directed toward defining proportico of Bound 
fioldo which oecn to have a perceptible and quantitative 
effect upcti the Goodneso of the acoustical environment. 

Qno ouch proporty io ^diffusion.* Without delving into 
tho technical, psycho-acoustic aepecto of tho problem, 
it vrculd be advisable to diocuso same of tho more intui- 
tively perceptible aspects of acoustics which aro 
attendant to the property of diffusion. 

A sound field is defined as being completely diffuse 
if it has uniform energy density within tho region con- 
sidered, and if the direction of proposition of tho wavo 
fronts at any arbitrary selected points are wholly rand cm 
in distribution. In discussing came of tho trends in the 
study of sound waves In rooms^" Morse and Bolt say, "It 
now appears that reverberation tine alone io not always a 
sufficient noasuro of auditorium excellence. It is 
desirable to have a moan square pressure as nearly uniform 
os posolble over the coating area. It io alco important 
to have at least a certain percentage of the sound reach 
the listener directly from tho speaker, and loss then a 
certain percentage roach tho listeners indirectly, after 
reflection from any single surface; ..." In tho light 



of the definition of dlffunenesc, the above quotation 
proscribes that for any sound field, thoro ia an optimum 
docree of diffusion \diich. produces a good acoustical en- 
vironment. Che may reason that .thla optimum degree of 
diffusion vd.ll depend upon the uno to which the Bound field 
io tc bo adopted i.o. lecture hall, broadcast otudio, 
cinema, music hall* and more particularly, upon the 
acoustical tradition of tho population in the Bound field. 
For instance, tho population of Britain may bo said to have 
slightly different preferences than the population of tho 
United States. Those preferences would probably bo based 
upon tho way the British have been used to hearing their 
music or lectures and also upon tho peculiarities of the 
language. Tho degree of diffusion in any sound field will 
in general bo a function of tho location of tho source, the 
shape of the boundaries enclosing the sound fiold, and the 
shape, absorption coefficient and distribution of the sur- 
faces in tho sound field. 

Ac a sort of intuitive example of tho effects of 
diffusion lot us visualise a full symphony orchestra play- 
ing in an open fiold. Generally, the conductor would be 
able to hoar clearly only those instruments which wore 
near to him or whose radiation was directive and boomed 
in his direction. Tho musicians would bo ablo to hoar 
proportionally less of the total effect being produced by 
the orchestra as a group. 



Under theoo circumstances, it io probable that a 
very poor porfornonco would bo tho result of this lack 
of dlf fusion or nixing. On the other hand, if cay a 
email c trine onsomblo was performing in a very diffuse 
broadcast studio, oach aombor of tho ensemble would bo 
able to hear tho combined of feet of tho Group, but ho 
would hoar relatively little of tho sound that hie par- 
ticular instrument was producing. Such a condition would 
again result in an inferior porfornanco. With tho proper 
degree of diffusion, however, it io thoorisod that con- 
ditions would bo optimum for tho production of tho desired 
acoustical offoct. 

In a technical report Number B 058, serial No. 1953/ 
29, dated October 1953> tho British Broadcasting Company 
discussed tho property of diffusion of sound fiolds. It 
was said that imperfect diffusion nay appear in at loaot 
two measurable forms, (1) as a change in pros sure with a 
change in position, tho frequency of tho sound remaining 
constant; (2) as a change in pressure with changoo in fre- 
quency, tho positions of source and receiver remaining 
unchanged. It was further stated that Imporfoet diffusion 
Indlcatod in tho steady state measurement by variation of 
pros euro with frequency or position is shown by deviation 
in tho pressure, time relations from an exponential decay 
curve. These statements oro perfectly reasonable; however, 
tho functional relationship between tho degree of diffusion 
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and the tranciont response of a sound field Is not ap- 
parent* 

Correlation techniques seen to offer a method of 
obtaining a moaauro of tho diffuoion of a sound field 
which do pond o diroetly upon tho definition of diffuoion. 
Moreover, this method employ a almost the camo equipment 
and arrangements as are used for tho transient ro spans© 
raoaourenente . 

4.2 COKRK.ATXON METHOD OF DIFFUSION MEASUREMENT 

Lot uo visualize tho following expo rlnont . 

V© shall placo a microphone in the sound field created 
by a loudspeaker driven by a wide band noise voltage 
generator. 

Nov; let us auto-eorrelato the output voltage of tho 
microphone. If the noise voltage has an Infinitely wide 
band flat spectrum, then it has boon shown? that tho auto- 
correlation function of tho microphone voltage is cscactly 
tho combined impulse reopens© of the loudspeaker end the 
microphone. For noloo voltages of finite bandwidth, this 
auto- correlation function is tho transient response of tho 
system to a pulse having tho corresponding spoctral cam- 
position. In general, this auto-correlation function has 
a maximum peak for zero time delay end a first axis cros- 
sing at a time delay dopondont upon tho bandwidth of tho 



noino voltage oo - Qx and tho quantity oo + mi. A 
aorlen of ouch auto-correlations Is ohov/n in Fig. 4.1. 

As explained In Chapter I, tho bandwidth and the 
first axis crossing of tho auto-correlation function vary 
inversely. That is, tho wider tho bandwidth, tho chortor 
tho difforence in delay times bet’./eon tho poak of the cor- 
relation function and tho first axis crossing. It is 
possible to measure the time delay corresponding to this 
first axis crossing with a high dogroo of accuracy. Tho 
tino delay register associated with tho analog correla- 
tion conputor poinits readings of tiao dolay to tho tenth 
of a millisecond. It is possible to interpolate botwoon 
those roglster readings by tho following method. 

a) Into oach channel of the time delay, food a 
standard 10 hops signal voltage. 

b) Measure tho phase dlfforence of tho outputs of 
tho tino delay for this 10 keps signal. 

c) Since the period of the 10 keps signal is pre- 
cisely 1/10 of a millisecond, it is possible to 
hand adjust tho time dolay mechanism to a value 
corresponding to any fraction of a period in 
torso of degrees of phase difference between the 
tiao delay outputs. 

d) Tith the time delay mechanism thus adjusted, one 
can now remove the 10 keps signal, apply the noise 
voltage, and measure the voltage output correspond- 
ing to tho value of the auto-corrolatlon curve for 
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that tlno co lay. 

How, knowing the time delay of the axis crossing, 
one can place two such identical microphones in tho 
sound field of a source radiating spherical sound waves. 
If the microphones are placed a distance apart equal to 
the time delay to the firGt axis crossing of the auto- 
correlation curve, and then orientated in such a way that 
tho sound fronts from tho source strike both nikos at the 
sane instant (Fig. 4.2 ) ; then tho cross- correlation 
function at zero tlno delay for the output voltage of 
the two microphones will correspond to tho value of the 
auto-correlation curve for = 0. Howovor, if the source 
microphone array is turned 90 degrees in a horizontal 
plane co that the sound waves strike one microphone, is 
subjoctod to a physical delay corresponding to tho axis 
crossing tino, and then strikes the second microphone 
(Fig. 4.3), thGn the corresponding cross-correlation 
function will be aero. Thus it is possible to obtain a 
plot of cross-correlation vs the orientation the micro- 
phones array in degrees from 0 to 90°. If wo cpccify 
that zero degrees corresponds to the orientation when 
both mikes are equidistant from the source, then to can 
convert our previously measured accurate plot of the 
auto- correlation curve vs tine delay to a plot of cross - 
correlation function vs orientation in degrees by means 
of the formula = d sin 9 where d is the niko 
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FIs* 4.2 Hiotosraph of microphone array oriented 

at 0 degrees in oound field of loudspoa&or. 
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Pig. 4,3 



Riotograph. of microphone array oriented 
at 90 degrees in sound field of loudspeaker. 
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separation in feet, c is the Bpeed of sound in feet per 
millisecond, and is time delay in milliseconds. These 
calculations have been plotted for a noise spectrum flat from 
100 cps to 10 kcps, having a roll off of 6 db per octave 
beyond those frequencies (Pig. 4.4). 

The microphone array together with the multiplying 
and integrating circuits of the analog correlator comprise 
an instrument which gives a measure of the properties in- 
volved In the definition of diffusion. 

Suppose the microphone array was placed in a perfectly 
diffuse field. Then, no matter which way the array was 
turned, the value of the cross-correlation function would 
be the same, since by definition, a diffuse field is one 
in which the wave fronts appear to be coming uniformly 
from all directions . It is clear that for variations in 
diffusion, the related cross-correlation curve would have 
some characteristic shape. 

4.3 EXPERIMENTAL RESULTS 

EXPERIMENT 1 

Two Altec 21-BR-150A microphones were placed two 
inches apart in arechoic space in the sound field pro- 
duced by an Altec 728-B closed -box -baffled loudspeaker. 

The cross correlation of the clipped output voltage of 
the microphones wa3 plotted as a function of orientation, 

Pig. 4.5, from 0 degrees. Pig. 4.2, to 90 degrees, 
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Fig. 3. Tho curve in Fig. 4.5 13 a measure of the dif- 
fusion of anacholc space. As should be expected, the 
source shows that the sound field produced by the loud- 
speaker in anechoic space is not diffuse . 

The valuo of microphone separation used in this first 
experiment was obtained by trial and error. It is impos- 
sible to got tho desired results by spacing the micro- 
phones the theoretically correct distance apart because 
of the finite size of the microphone diaphragms. 

EXPERIMENT 2 

The following experimental work was performed in the 
model studio, Room 20-F-009A, located in the Acoustics 
Laboratory at M.I.T. This studio was designed as a 
diffuse room. A physical description of this room is 
contained in reference 12. 

With the loudspeaker facing the wall shown in Figs. 
4.2 and 4.3> diffusion data was taken with the microphone 
array placed at various distances for a line perpendicular 
to the radiating face of the loudspeaker enclosure. The 
noise voltage and other equipment was the same as in 
Experiment 1. As in Experiment 1, the microphone output 
voltagos were clipped 20 db below their rms value. 

The results of the diffusion measurements together 
with the plot of the transient response of the studio are 
shown in Figs. 4.6 and 4.7 respectively. 
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EXPERIMENT 3 

Tho following experimental data was taken in Room 
10-390-A, the hard plaster room. In these experiments 
tho data was tape recorded and then brought into the 
laboratory and reduced. 

Tho carao equipment used in Experiments 1 and 2 was 
used in the following tests. 

Diffusion data and transient response data was taken 
for various conditions of the room. For these toots, the 
loudspeaker was placed facing a comer of tho rocn and the 
diffusion data was taken for various positions in the room. 
The positions wore located at distances of 20, 15, 10, 8, 
and 5 feet from the loudspeaker on a line between tho 
loudspeaker location and tho diagonally opposite comer 
of the roam. The transient response data was taken at the 
15 foot position. 

The above described data was taken, reduced, and 
plotted for the following room conditions: ' 

a) All walls of the room were bare. Figs. 4.8 and 

4.9. 

b) One long vrall of the rocrn was covered with a 
highly absorbent Fiberglas filled quilt. Figs. 

4.10 and 4.11. 

c) One long vrall and one short vrall were covered with 
quilting. Figs. 4.12 and 4.13. 

d) Two short walls and one long vrall were covered with 
quilting. Figs. 4.14 and 4.15. 
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Figure 4.9 
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Figure 4.15 
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4.4 DISCUSSION OF EXPERIMENTAL RESULTS 

In reviewing the experimental results of Section 
4*3 it appe are that the measuring system and the tech- 
niques employed yield results similar to those obtainable 
with a highly directional microphone. It appears that 
the practice of setting the microphones a. certain dis- 
tance apart corresponds to making a microphone which is 
extremely sensitive in measuring the direction of prop- 
ogatlon for a sound having a particular spocular 
composition. Generally, setting the microphones further 
apart makes the array more sensitive to spectra having 
lower geometric mean frequencies . Insufficient data is 
available to permit an evaluation of this method of 
measuring diffusion. However , the author believes that 
the data presented herewith is sufficient to indicate a 
possible experimental method of obtaining a measurement 
of the degree of diffusion of a sound field. Moreover, 
it Is believed that this method warrants further study. 

4.5 SUGGESTIONS FOR ADDITIONAL WORK 

On the basis of the discussion and measurements of 
diffusion made in connection with this project, it is 
felt that there is some potential merit in the techniques 
employed. As a suggestion for additional work, it is pro- 
posed that some consideration be given to further diffusion 
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measurements made T&th various bandwidtho of nolao voltage. 
Such a project vould make extensive experiments in some 
room such as the model studio. One possible series of 
experiments would involve measurements for various micro- 
phone separations (selected in a manner similar to that 
of Experiment 1, Section 4.3) and various conditions of 
diffusion (adjustments in diffusion to be made by 
rearranging variable diffusers in the model studio) . 
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V. CONCLUSION 

5.1 SUMMARY OP CORRELATION MEASUREMENTS 

Correlation measurements as studied in this project 

t 

have been of two types - (1) measurements of the transient 
response of rooms and (2) measurements based on the 
definition of a diffuse sound field. It is an interest- 
ing coincidence that both types of measurements can be 
made with practically the sane experimental setup. 

Transient response data may be taken by cross correlating 
between the microphone output voltages of two microphones, 
one of Which is located near the source and the other in 
the reverberant field of the room. Diffusion measurements 
are mado in the same way, by cross correlating between the 
microphone output voltages of two microphones which are 
spaced a particular distance apart and oriented in a 
certain way with respect to the source. 

The above illustrated coincidence serves to point out 
the versatility of correlation techniques. Unfortunately, 
this coincidence also serves to illustrate one of a number 
of pitfalls that one may encounter in the use of such 
techniques. Improper microphone placement in either type 
of measurement would lead to spurious results. 

It was stated before that diffusion measurements by 
correlation techniques seemed to yield results similar to 
those obtainable by use of a highly directional microphone. 
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It should be of Interoot to point out that spacing 
microphones a distance apart which is calculated from 
auto- correlation data corresponds to calibrating the 
microphone array for maximum directionality for the 
particular frequency spectrum of noise voltage being 
used. 

This directive microphone array constructed by cor- 
relation techniques appears to havo two major lobes in 
its directivity pattern. These major lobes are 180 
degrees apart. The calibration process consists of 
spacing the microphones so as to eliminate any minor 
lobes in the directivity characteristics of the array. 

It should be pointed out that the directivity which 
we have been talking about refers only to the ability of 
the microphone array and the associated correlation equip- 
ment to discriminate between the presence of a plane wave 
front which roaches both microphones at the same instant 
and a plane wave front which reaches each microphone at a 
different instant. 

In a previously mentioned report by the Eritlch 
Broadcasting Company, it was concluded on the basic of 
their experimental studies that the methods of short pulse 
analysis are not satisfactory for the investigation of 
diffusion in full scale rooms. The pulse methods refer- 
red to in tills conclusion involved - (a) studies of the 
irregularity of the envelope of the decay for a short 
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pule© of sound and (b) variation in thoso short pule© 
decay irregularities with position in tho room. 

Another group of investigators J U3ing a direc- 
tional microphono technique attempted to obtain a Gteady 
state moaouro of tho diffusenoss of sound fields. The 
conclusion to their research cited tho necessity for a 
microphone that could be made highly directional over 
the whole audio frequency range. 

It appears that the correlation method of measuring 
transient response is suitable for use in full scale 
rooms. Unlike short pulse methods, correlation analysis 
of transient response may bo made to yield short term 
decay rates for various bandwidths of frequencies. On 
the basic of this latter fact, it may be concluded that 
correlation analysis presents a connecting Unit between 
room response data gathered by steady state methods and 
that gathered by pulse methods. 

In addition, it can be seen that correlation tech- 
niques, when used in conjunction with a microphone array 
calibrated for tho purpose, can be made to yield data 
similar to that obtained by directional microphone 
methods. Further, tho microphone array can bo calibrated 
for any bandwidth by use of data taken from its auto- 
correlation curve. 

From these facts it may be further concluded that the 
correlation method of analysis of bounded sound fields is 
potentially more informative than either steady state 
methods or pulse methods taken alone. 



